To be published in the Astrophysical Journal 



FAST ROTATION OF NEUTRON STARS 

Jean-Pierre LASOTA \ Pawel HAENSEL 2 , 
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Marek Artur ABRAMOWICZ 3 
ABSTRACT 

We show that for realistic equations of state of dense matter, the universal 
proportionality factor relating the maximum rotation rate of neutron stars due 
to mass-shedding limit to the mass and radius of maximum allowable mass 
configuration of non-rotating models results from a universal proportionality 
between masses and radii of static maximum-mass neutron stars and those of 
maximally rotating configurations. These empirical relations cannot be obtained 
in the slow rotation approximation. 

Subject headings: pulsars-stars: neutron-stars: rotation 



1. INTRODUCTION 

The question of how fast a pulsar can spin has recently been discussed by several 
authors (see e.g. Friedman & Ipser 1992, Cook et al. 1994a, b). Independently of the 
possible mechanisms of pulsar spin-up, for any particular equation of state there is always 
an upper limit for the final allowed rotational speed, fi max . The limit is determined either 
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by the mass shedding condition, i.e. that the angular velocity of the configuration, stable 
with respect to axisymmetric perturbations, equals the Keplerian velocity at the surface, or 
by the condition of the onset of non-axisymmetric (e.g. gravitational radiation reaction) 
instabilities. In this paper we discuss the limit given by the mass shedding condition. 

Haensel and Zdunik (1989; hereafter HZ) noticed that for realistic equations of state 
of dense matter, the numerically calculated values of the shedding limit fi max can be fitted, 
with an accuracy better than 5%, by an empirical formula 



where M max is the maximal mass of the non-rotating neutron stars with the same equation of 
state, -R max is the radius corresponding to M max , and Cq is a dimensionless phenomenological 
constant, independent of the equation of state. fi max is an angular velocity of rigid rotation 
as measured by a stationary observer at infinity. HZ determined that the best fit is for 



The calculations of Friedman, Ipser and Parker (1989), performed for a very broad 
set of equations of state, yielded Cq = 0.66, quoted in Friedman (1989), Friedman and 
Ipser (1992). It should be mentioned, that the value of Cq = 0.62, quoted in the original 
paper of Friedman, Ipser and Parker (1989), was actually a very rough estimate of Cq (J.L. 
Friedman, private communication). Calculations based on recent very accurate numerical 
methods are in good agreement with the original HZ choice of Cq. For example, values 
of O max calculated for several equations of state by Lattimer et al. (1990) differ from the 
HZ version of the empirical formula by less than 4%. Most recent calculations in the full 
framework of General Relativity by Salgado et al. (1994a,b; hereafter SBGH) based on the 
spectral methods (Bonazzola et al. 1993) show that configurations with causal equations 
of state (EOS) satisfy (1.1) with Cq = 0.67 with an accuracy of better than 5% (Haensel, 
Salgado and Bonazzola 1995, hereafter HSB). Results obtained for causal EOS by Cook 
et al. (1994b) lead to very similar "best fit" value of Cq. HSB found that the empirical 
relation with Cq = 0.67 fails for configurations constructed with non-causal equations of 

1 /2 

state, where the sound speed, {dP/dp)g , may be greater than the speed of light within the 
neutron star models. 

Empirical relations of universal character, valid for a broad range of realistic equations 
of state of dense matter, can obtained also for other bulk parameters of neutron star 
models. For example, a simple and surprisingly good universal relation, connecting the 
maximum moment of inertia for slow, rigid rotation, J max , to the mass and radius of a static 
configuration with maximum allowable mass (which is different from that with J max !) was 
pointed out by Haensel (1992). The existence of such universal empirical relations is of 




(1.1) 



Cq = 0.67. 
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practical importance: for any realistic EOS, empirical formulae enable rapid and still quite 
precise calculation of the upper bounds on global parameters of neutron stars from the 
easily calculated parameters of the static configuration with maximum allowable mass. A 
different type of a universal formula was obtained by Ravenhall and Pethick (1994). Their 
formula, valid for a broad range of realistic equations of state of dense matter, and useful 
for all except lightest neutron stars, expresses the moment of inertia in terms of stellar mass 
and radius. 

The few attempts to explain the empirical relation, Eq. (1.1), were not concluded 
with satisfactory results. The most elaborate discussion was presented by Weber and 
Glendenning (1991; 1992). They used numerical models of rotating neutron stars calculated 
in the slow rotation approximation of Hartle and Thorne (1968) to show that these also 
obey the empirical formula, albeit with Cq ~ 0.75. Although this is obviously an important 
result, still lacking is a clear physical understanding of the problem. 

In a recent work Glendenning and Weber (1994) derived a formula which relates ^ max 
to M max /i?^ ax , in the slow rotation approximation, where M max and -R max correspond to 
rotating configurations, but they do not discuss the connection between this relation and 
the formula (1.1). As we shall see below the slow rotation approximation cannot be applied 
consistently to maximally rotating neutron stars, as pointed out already by Hartle and 
Thorne (1968). 



2. THE EMPIRICAL RELATIONS BETWEEN NON-ROTATING AND 
MAXIMALLY ROTATING CONFIGURATIONS 



One can see that, to a good approximation, the "constant" Cn is in fact a slowly 
varying function of one variable, characterizing the EOS of dense matter, C^(EOS) (Table 
1 and HSB). Indeed Fig. 1 shows the values of Cq as given by eq. (1.1) for 12 maximally 
rotating models calculated by SBGH, as a function of the parameter 

2GM max . . 

■''max 1 ' 

for static maximum mass configurations. It is clear that C^(EOS) is slowly varying, and for 
the range of parameters of interest, can be well approximated by a monotonic function of x s . 
If one restricts to realistic EOS, which are both causal and stiff enough to support observed 
masses of pulsars, the range of relevant x s becomes rather narrow and the approximation of 
Cn(EOS) by a constant is quite satisfactory (see HSB). 
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It is well known that for realistic EOS, rotation increases the value of maximum mass 
of a neutron star by about 20%. We find however that, with a very good approximation, 
there exists a universal relation between the maximum mass of non-rotating neutron star 
configuration and the mass of a configuration rotating with an angular speed f2 max - The 
relation is: 



where M max (rot) and M max (stat) are respectively the mass of the configuration in maximum 
rotation and the maximum mass of the static neutron star for the same equation of state. 
Of course, a similar relation holds trivially for any specific EOS, and yields a specific value 
of the proportionality constant, Cm(EOS). The values of Cm(EOS), calculated for a broad 
set of EOS using the results of SBGH, are given in Table 1. The best fit of relation (2.2) to 
numerical results of SBGH is obtained for C M = 1-18. Then, relation (2.2) reproduces exact 
results within better than 3% (Fig. 2). 

A similar relation is found for the radii of the corresponding configurations: 



where _R max (rot) and _R max (stat) are the radii of the maximally rotating and the static 
configurations respectively. The best fit value of a constant Cr in relation (2.3), based on 
numerical results of SBGH for causal EOS, is Cr = 1.34. Relation (2.3) reproduces then 
exact results for causal EOS within better than 4% (Fig. 3). 

As in the case of relation (2.2), one can introduce the factor Cr(EOS) (Table 1). The 
dependence of Cr on the EOS can be well approximated by a monotonically decreasing 
function of x s - this is visualized in Fig. 4. No such a trend was found for C M (EOS)), Fig. 
5. 

Our analysis have been based on a set of numerical results, obtained for twelve realistic 
EOS in SBGH. We have restricted to the SBGH set in order to keep the homogeneity of 
the sample of numerical results. Numerical results obtained by other authors (Friedman, 
Ipser & Parker 1989, Cook et al. 1994b) are found to be slightly different from those of 
SBGH. This is due to the fact, that the precise determination of the maximally rotating 
configuration - which requires very high precision of numerical procedure - turns out to 
be sensitive to such details as, e.g., the interpolation method used for the determination 
of the maximum frequency model (Stergioulas & Friedman 1995, E. Gourgoulhon, private 
communication). However, uncertainties resulting from this dependence on the specific 
sample of numerical results obtained for realistic EOS, are consistent with the precision of 
relations (2.2), (2.3), discussed above. 



M max (rot) = C M M max (stat) 



(2.2) 



#max(rot) = C Rj R max (stat), 



(2.3) 
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One can easily see that the empirical constant may be obtained, within a very good 
approximation, from the formula: 

As can be seen from Table 1, for causal EOS Eq. (2.4) reproduces the actual values of 
C^(EOS) within better than 5%. In the case when one considers also non-causal EOS, the 
precision of approximation (2.4) worsens to 7%. 

Relation (2.4) implies, that the angular velocity for mass shedding is approximated by 
the formula 

The value of f2 max can be thus well approximated by the frequency of a particle in stable 
circular orbit at the equator of a fictitious non-rotating star of mass equal to M max (rot) 
and of the radius equal to the equatorial radius of the maximally rotating configuration, 
i? max (rot) (in the case of non-rotating star the general relativistic formula for the particle 
frequency is identical to the newtonian one, see e.g. Misner, Thorne & Wheeler 1973). 

Up to now, we considered only realistic EOS of dense matter. It is instructive to study 
also the case of polytropic configurations with an equation of state in the form P = Kn^ r , 
where is baryon density of matter. Extensive calculations of rotating polytropic models 
were presented in Cook et al. (1992, 1994a). For a fixed value of T, static and rotating 
models exhibit useful scaling properties with respect to change of K (Cook et al. 1992, 
1994a). Consequently, the relativistic parameter x s , and Cm(EOS), Cr(EOS) turn out to be 
functions of T only; they are given in Table 2. For a fixed value of T, relations (1.1), (2.2), 
and (2.3) are thus exact, with T-dependent values of the numerical coefficients. 

The realistic EOS are not polytropes, and their local adiabatic index, 
r = (n h /P)dP/dn h , depends on the density, T = T(n h ). Two specific examples of 
the density dependence of T are shown in Fig. 6. In both cases, the value of T varies within 
the relevant interval of nb by more than 30% of its maximum value. Clearly, such EOS 
cannot be represented by single polytropes. This explains, e.g., the lack of monotonic trend 
in the dependence of Cm(EOS) on x s (EOS) for realistic EOS, displayed in Fig. 5. Such a 
non-monotonic, irregular behavior, characteristic of realistic EOS, is to be contrasted with 
a monotonic dependence Cm(x s ) for polytropes, Table 2. In spite of this irregular behavior 
of Cm(EOS) for realistic EOS, the values of C(EOS) for realistic EOS, calculated from (2.4), 
show a clear trend for a monotonic increase with x R . This results from a clear trend in 
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Cr(EOS) to decrease with x s , which - magnified by the third power within the bracket of 
formula (2.4) - dominates the C - x s relation. 

The trend for a monotonic increase of Cn with x s , Fig. 1, can be well reproduced by a 
linear function C l Q n (x s ) = 0.468 + 0.378:r s . This function, obtained by a least squares fit to 
the points shown in Fig. 1, reproduces C^(EOS) (and f2 max (EOS), if inserted in Eq. (1.1)) 
with precision better than 1.5%, with typical relative error being less than 1%. It should 
be stressed, however, that such a linear approximation C^ n (:r s ) is valid only within a rather 
narrow interval of x s , characteristic of realistic EOS, 0.45 < x s < 0.7. In contrast to the 
one-parameter empirical formula with Cq = 0.67, C^ n (x s ) does not work for the EOS of the 
free neutron gas (see HSB). 



3. THE SLOW ROTATION APPROXIMATION 



Weber and Glendenning (1991, 1992) tried to derive the empirical formula, equation 
(1.1), by using the slow rotation approximation of the Einstein equations. One should 
notice however that a neutron star rotating with fl max may not be really considered to be a 
slow rotator (Hartle & Thorne 1978). Indeed, let us define a dimensionless angular velocity 

«.'-"V(^). (3-D 

where M is the mass of the static star, R its radius and Q is the angular velocity measured 
by observers at infinity. The assumption of the slow rotation means that (Hartle 1967), 

« i. (3.2) 

From equation (1.1) with Cq = 0.67 it follows that 

^ 2 max « 0.45, (3.3) 

so that the use of slow-rotation approximation for maximally rotating configurations 
requires some additional justification. It is unlikely that this approximation will give 
accurate results considering the fact that at the stellar surface the linear speed of rotation 
is a significant fraction of the speed of light (Hartle & Thorne 1968): 



v s i f x s y/ 2 



c \/2\l-x t 



(3.4) 
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For the maximally rotating models of SBGH one gets typically v s /c « 0.7. In the slow 
rotation approximation, one can obtain for f2 max an expression of the form, which seems to 
be similar to that of the empirical formula (1.1) (see also Glendenning and Weber 1994): 

Onax = C ( ^ r x(rot) \ 1/2 + 0{n% (3 _ 5) 



with 



r 



#max( r0t ) 



1 R G o , 1 { R \ A {, R G\ '2/..N {-, Q M ° 2 



Mi? 2 i? v ~ MR2 \ r g/ V i? / zy ' v J2 

Here Rq is the gravitational radius, 



(3.6) 



Ro - ™ (3.7) 



C 2 



J is the angular momentum, I the moment of inertia, Q is the quadrupole moment of the 
rotating configuration, and Q'^u) is the derivative, with respect to u — 1 — R/Rg, of the 
associated Legendre function of the second kind. The quantities M and R are those for the 
maximally rotating configuration, but - within our approximation — we can as well replace 
them by M max (stat), i? max (stat). 

Equations (3.5) and (3.6) give f2 max in terms of the mass, equatorial radius, moment 
of inertia, angular momentum and the quadrupole moment of the maximally rotating 
configuration. At first glance (and neglecting terms ~ f^ 3 and higher), expression (3.5 - 
3.6) may seem to be similar to empirical formula (1.1). It has been shown by Abramowicz 
and Wagoner (1978) and recently confirmed by Ravenhall and Pethick (1994) that moments 
of inertia, expressed in the units of MR 2 , are - to a good approximation - functions of 
only x = R/Rg- Moreover, results of SBGH show that QMc 2 / J 2 is a decreasing function 
of x s . Expression (3.5) seems thus to possess similar properties as the empirical formula. 
Unfortunately it is not the empirical formula since it involves M max and R max of rotating 
configurations, and therefore it corresponds to some approximation of eq. (2.5) and not to 
eq. (1.1). Further expansions of M max (rot) and _R max (rot) in are not justified - in view 
of the large value of this parameter (notice that rotation increases the value of -R max by 
some 30%). So, only some external "empirical input" (such as assumption of a "typical" 
effect of rotation on M max , -R max , and on eccentricity of rotating star, made in Weber & 
Glendenning (1992)) can lead to an expression of type (1.1). A consistent application of the 
slow rotation approximation cannot reproduce the empirical formula for fi max . This is not 
very surprising if one considers that one should get to the 8-th order in to get a precision 
of 4%, characteristic of empirical formula. 
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4. DISCUSSION AND CONCLUSIONS 

The proportionality constant appearing in the empirical formula for Q max for realistic 
EOS of dense matter, is in fact a function of the relativistic parameter x s . In the range 
of parameters describing maximally rotating configurations with causal equations of state 
Cq is a rather slowly varying function of x s , which results in a high precision of empirical 
formula with an appropriate choice of a universal proportionality constant. We found 
universal relations connecting maximal masses and corresponding radii of static neutron 
stars, calculated for realistic EOS of dense matter, with those of maximally rotating 
configurations. The empirical formula follows from those relations. 

Although the slow rotation approximation allows one to reproduce some of the 
properties of the empirical formula this approximation is not appropriate for maximally 
rotating configurations. 
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FIGURE CAPTIONS 

Fig. 1. Proportionality constant C^(EOS) for twelve realistic EOS versus relativistic 
parameter x s . Filled circles for causal EOS, open circles for EOS which are non-causal 
within massive neutron star models. 

Fig. 2. Mass of maximally rotating neutron stars versus maximum mass of non-rotating 
models for realistic EOS. Filled circles for causal EOS, open circles for EOS which are 
non-causal within massive neutron star models. Straight line gives the best fit to exact 
results obtained for causal EOS. 

Fig. 3. Equatorial radius of maximally rotating neutron stars versus the radius of the 
static maximum mass configurations for realistic EOS. Filled circles for causal EOS, open 
circles for EOS which are non-causal within massive neutron star models. Straight line 
gives the best fit to exact results obtained for causal EOS. 

Fig. 4. Proportionality constant Cr(EOS) for twelve realistic EOS versus relativistic 
parameter x s . Filled circles for causal EOS, open circles for EOS which are non-causal 
within massive neutron star models. Dashed horizontal line corresponds to the best fit to 
exact results for causal EOS. 

Fig. 5. Proportionality constant Ca/(EOS) for twelve realistic EOS versus relativistic 
parameter x s . Notation as in Fig. 4. 

Fig. 6. Adiabatic index, T, versus baryon density, rib, for AV14 + UVII (solid line) and 
UV14 + TNI (dashed line) EOS. The calculation of T has been based on the polynomial fit 
of Kutschera & Kotlorz (1993) to the tabulated EOS of Wiringa et al. (1988). 
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TABLE 1 

REALISTIC EQUATIONS OF STATE 



Equation of state 


x s 




Cm 


Cr 


C 


Glendenmng 1985 case 2" 


0.467 


0.648 


1.1741 


1.3960 


0.657 


/"I 1 1 * infir il -111 

Glendenmng 1985 case 1 


0.480 


0.659 


1.1783 


1.3794 


0.670 


Glendenmng 1985 "case 3" 


0.516 


0.658 


1.1746 


1.3740 


0.673 


Diaz Alonso 1985 model 11 


0.524 


0.666 


1.1704 


1.3575 


0.684 


T T T 1 ] "1 -I /~\ /~\ -1 1 (III T T "T T 

Weber et al. 1991 A£ u onn + HV 


0.533 


0.661 


1.1983 


1.3786 


0.676 


Bethe & Johnson 1974 model IH 


0.554 


0.672 


1.1597 


1.3409 


0.693 


Wiringa et al. 1988 model UV14 + TNI 


0.569 


0.687 


1.1826 


1.3277 


0.711 


Pandharipande 1970 


0.577 


0.684 


1.1633 


1.3226 


0.709 


Haensel et al. 1980 "model 0.17" 


0.614 


0.704 


1.2138 


1.3103 


0.7345 


Friedman & Pandharipande 1981° 


0.619 


0.701 


1.1807 


1.3145 


0.721 


Wiringa et al. 1988 model UV14 + UVII a 


0.656 


0.719 


1.1824 


1.2690 


0.761 


Wiringa et al. 1988 model AV14 + UVII a 


0.667 


0.722 


1.1904 


1.2671 


0.765 



a Non-causal within central cores of massive neutron stars 
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TABLE 2 

POLYTROPIC EQUATIONS OF STATE 
I x s Co. Cm Cr C 



2.00 0.438 0.624 

2.25 0.509 0.640 

2.50 0.563 0.660 

2.75 0.605 0.6745 

3.00 0.637 0.6923 



1.1504 1.4242 0.631 

1.1766 1.4032 0.652 

1.1961 1.3745 0.679 

1.2111 1.3540 0.6985 

1.2225 1.3288 0.723 



